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Abstract 

Four classic criteria used to the classification of complex flows are discussed here. These criteria are useful to 
identify regions of the flow related to shear, elongation or rigid-body motion. These usual criteria, namely Q, A, 
A 2 and Acr/Acz, use the fluid’s rate-of-rotation tensor, which is known to vary with respect to a reference frame. 
The advantages of using objective (invariant with respect to a general transformation on the reference frame) 
criteria are discussed in the present work. In this connection, we construct versions of classic criteria replacing 
standard vorticity, a non-objective quantity, by effective vorticity, a rate of rotation tensor with respect to the 
angular velocity of the eigenvectors of the strain rate tensor. The classic criteria and their corresponding objective 
versions are applied to classify two complex flows: the transient ABC flow and the flow through the abrupt 4:1 
contraction. It is shown that the objective versions of the criteria provide richer information on the kinematics 
of the flow. To cite this article: R. S. Martins, A. S. Pereira, G. Mompean, L. Thais, R. L. Thompson, C. R. 
Mecanique xxx (2015). 

Resume 

Une perspective objective pour des criteres classiques de classification d’ecoulements. Quatre criteres 
classiques utilises pour la classification des ecoulements complexes sont etudies ici. Ces criteres sont utiles pour 
identifier les regions de I’ecoulement liees au cisaillement, a I’extention ou au mouvement de corps rigide. Ces 
criteres habituels, a savoir Q, A, A 2 and Acr/Acz, utilisent le tenseur taux de rotation du fluide qui est connue 
pour varier par rapport au systeme de reference. Les avantages d’utiliser des criteres objectives (invariants par 
rapport a une transformation generale pour un systeme de reference) sont discutees dans le present travail. A 
cet egard, nous construisons des versions des criteres classiques en remplagant la vorticite standard, une quantite 
non-objective, par la vorticite effective, un taux de rotation par rapport a la vitesse angulaire des vecteurs propres 
du tenseur taux de deformation. Les criteres classiques et leurs versions objectives correspondantes sont appliques 
pour classifier deux ecoulements complexes : I’ecoulement ABC transitoire et I’ecoulement a travers une contraction 
brusque 4 :1. Les version objectives des ces criteres fournicent des informations plus riches pour la cinematique 
de I’ecoulement. Pour citer cet article : R. S. Martins, A. S. Pereira, G. Mompean, L. Thais, R. L. Thompson, 
G. R. Mecanique xxx (2015). 
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1. Introduction 

In Fluid Mechanics, flow visualization is an important subject, since fundamental aspects of the flow 
can be captured by observation. Post-processing Computational Fluid Dynamics (CFD) data is also a field 
that makes important contributions for the understanding of the flow. Complex flows exhibit different 
kinds of motion that depend on position and time. In these flows, it is common to find swirling motions 
in different parts of the domain. In order to locate and visualize these regions, a criterion of vortex 
identification is generally used to see the manifestation of the rotation character of the flow. 

However, the concept of a vortex is still cause for dissension within the scientific community. As a 
consequence, there a several criteria available in the literature that are used to identify rotation structures 
in the flow. In other words, there is no quantity, in the mathematical sense, that is consensually accepted 
in the literature as a definition for a vortex. Some of the non-consensual issues that are present in this 
context are if the vortex is an Eulerian or a Lagrangian entity and if it should be defined in a kinematic 
or in a dynamical basis. 

Comparisons among the different criteria are still a subject of investigation (e.g. [1]). An important 
point raised by Haller [2] is the requirement that a vortex should be an Eulerian invariant entity, i.e. 
invariant under arbitrary changes of the reference frame. This requirement affects the vortex concept, 
since, before that work, only the Galilean invariance was invoked to define a vortex [3]. We can stress here 
that the classic vortex definitions, such as the Q-criterion by Hunt et al. [4], the A-criterion by Chong 
et al. [5], the A 2 -criterion by Jeong and Hussain [3] and the Acr/Aci-criterion by Chakraborty et al [6], 
enjoy only Galilean invariance (i.e. they are invariant to constant velocity translating frames). 

The arguments to adopt an objective criterion for vortex identification are the following. Eirst, if one 
observer identifies a certain region as being a vortex while, for another one, this region is not a vortex, 
there is no reason to privilege the verdict stated by one observer with respect to the other. Secondly, we 
have to have in mind the advantages of building a criterion for vortex identification. One clear purpose 
of identifying a region as being a vortex is to connect the rotational character of the flow with another 
phenomenon besides the flow itself. It is consensual that processes like: the convection in a heat transfer 
problem, the degree of mixture of different fluids, the percentage of components due chemical reaction 
in a flow, the intensity of polymer stretching due to the flow, and other transport phenomena problems, 
cannot be observer-dependent. Hence, if a vortex is non-objective, the logic of cause-effect that could link 
the flow character with one of these measurers of the intensity that a certain phenomenon is occurring is 
weaker, when compared to an objective criterion. 

In the present work we employ objective versions of four classic criteria largely used in the literature. 
The classic criteria and their respective objective versions are analysed and applied for two benchmark 
cases, the transient Arnold-Beltrami-Childress (ABC) flow [14,15,16] and the flow through a 
4:1 contraction. 


Email addresses: ramon.martinsOpolytech-lille.fr (R. S. Martins), anselmo.pereiraOpolytech-lille.fr (A. S. 
Pereira), gilmar.mompeanOpolytech-lille.fr (G. Mompean), laurent.thaisOpolytech-lille.fr (L. Thais), 
rthompsonOid.uf f .br (R. L. Thompson). 
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2. Classic criteria 


In the following we briefly present four criteria that are currently used in the literature to classify 
different regions of the flow. These criteria are Eulerian and Galilean-invariant and were recently selected 
by Pierce et al. [1] to evaluate for instance boundary layer flows. 

The Q-criterion was proposed by Hunt et al. [4] in the context of incompressible flows. Besides local 
pressure minima, they required that, to identify a vortex, the second invariant of the velocity gradient 
tensor, Vn (defined by Vn = {dujldxi)eiej^ where u = UiCi is the velocity vector field), should be 
positive. This condition can be expressed for incompressible flows as a function of the Euclidean norms ^ 
of the symmetric, D = (T + T^)/2, and skew symmetric, W = (L —L^)/2, parts of the velocity gradient 
(where L = {duildxj)eiej is the transpose of Vn). The condition for the Q-criterion can be expressed 
as follows. 


Q = \{\\Wf-\\Df)>0. (1) 

The A-criterion proposed by Chong et al. [5] is based on the assumption of an equivalence between 
a vortex and complex eigenvalues of the velocity gradient tensor^. Complex eigenvalues of the velocity 
gradient is a sign of vorticity dominance with respect to rate-of-strain, since the symmetric rate-of-strain 
tensor can only have real eigenvalues. Mathematically, the A-criterion can be defined as 


A = 



V 2 


2 

> 0 , 


( 2 ) 


where det(-) is the third invariant (determinant) of a given second order tensor. 

The A 2 -criterion proposed by Jeong and Hussain [3] is based on the idea of joining the local pressure 
minima condition to a vorticity predominance over the rate-of-strain in the same mathematical condition. 
By making some assumptions, neglecting some terms, this condition leads to 


As = Af < 0 , (3) 

t^2 I ti^2 o 2 

where A 2 ^ is the intermediate eigenvalue of the tensor D -\-W . 

The Acr/Aci-criterion proposed by Chakraborty et al. [6] is based on the concept that material points 
that follow orbits which remain compact during the revolutions around each other are in a vortex. Xd and 
Acr are the imaginary and real values, respectively, of the conjugate complex eigenvalues of the velocity 
gradient: 


\ \D+W \L 

^cr _ ^cr _ ^cr /a\ 

A. ■ A^+^ - ^ ^ 

The criterion was proposed in the form Xcr/^d < ^ and, therefore, a threshold parameter (e) is needed. 
This condition can only be applied in regions where there are complex eigenvalues of the velocity gradient. 
This criterion was conceived to be applied together with the Xd > S criterion proposed by Zhou et al. [7], 
where another threshold is needed. Since we will use a normalized dimensionless version of the criterion, 
these two conditions can be condensed into a single one, where to assume that Xd is significant is equivalent 
than to assume that Xcr/Xd is small. 


1. The Euclidean norm of a generic second order tensor A is ||A|| = ^tr(A • A^), where tr(-) is the first invariant 
(trace) of a given second order tensor. 

2. The expression velocity gradient is used interchangeably for L or Vit. 
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2.1. Objective redefinition for classic criteria 


Drouot [8], Drouot and Lucius [9] have shown that the relative rate-of-rotation tensor, W^ defined by 

W = w-n^, (5) 

is an objective quantity. In Eq. (5), accounts for the rate-of-rotation of the eigenvectors of D, and is 
defined by 


n 


^ = e?'e? 


(6) 


where eP are the normalised eigenvectors of and eP are their time derivatives. This quantity was 
used by Astarita [10], to form an index defined by the ratio of the norm of W to the norm of D. This 
index was used as means for flow classification dividing the domain into extension-like motions and rigid- 
body-rotation-like ones. Astarita [10] was seeking for an objective quantity when he proposed that flow 
classifier. 

Although a different path was followed by Dresselhaus and Tabor [11] and Tabor and Klapper [12] 
investigating the alignment and stretching of material filaments^ in an approach of dynamic systems, 
they have also arrived on the necessity of expressing these quantities with the help of the relative-rate-of- 
rotation tensor, which was called effective vorticity. The physical interpretation of the effective vorticity 
is discussed next. 

The more common interpretation of the vorticity tensor is that a vorticity component associated to a 
certain plane is the (arithmetic) mean of the rate-of-rotation of two filaments initially orthogonal to each 
other in that plane. However, a less adopted interpretation is that the vorticity is the rate-of-rotation of 
the filaments which are aligned to the eigenvectors of D. This fact induces the necessity of evaluating the 
vorticity with respect to the rate-of-rotation of the eigenvectors of D. When the eigenvectors of D are 
fixed in certain frame of reference, the vorticity tensor is responsible for a deviation between the filaments 
that are initially aligned with the eigenvectors and their respective directions. However, when the flow is 
complex and the eigendirections of D do change in time, the effective vorticity is the entity responsible for 
this deviation and, by consequence, for exposing different material filaments to the eigenvectors. This also 
means that when effective vorticity vanishes, the same material filament is exposed to the eigendirections 
of the rate-of-strain tensor. In this case, the filament aligned to the eigenvector corresponding to the 
highest positive eigenvalue is persistently stretched. As shown by Dresselhaus and Tabor [11], effective 
vorticity plays an important role on the dynamics of material lines and vorticity lines, affecting vortex 
stretching also. 

As suggested by Haller [2], an alternative but analogous form of the Q-criterion can be built by replacing 
vorticity by effective vorticity in Eq. (1) leading to 


0 = -II-DII") >0 


( 7 ) 


A similar procedure can be adopted for adapting the other criteria into an objective backbone. Therefore 
we can define the new versions of A-, A 2 -, and Acr/Aci-criteria, respectively, as 

^ fdet{D + W)\^ 



V 


V = Ar+^ 


■ 


> 0 


< 0 


(8) 

( 9 ) 


3. Material filaments being straight lines of infinitesimal size in the fiuid which can rotate, stretch, compress, but not 
bend. 
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( 10 ) 


Acr ^ Ag+^ ^ ^ 

\ , “ ^D+W ~ \L 

\i ^ci 

We highlight that all objective entities hereafter will be displayed with a hat. 


3. Results 


This section presents the results obtained by applying both objective and non-objective criteria to two 
cases: (1) the three-dimensional analytical flow field known as the ABC flow in its unsteady version, and 
(2) an abrupt 4:1 planar contraction. 

It is worth noting that all criteria have been normalized in order to produce values between 0 and 
1. Moreover, normalized values greater than or equal to 0 and less than 0.5 represent swirling-like or 
elliptical regions, whereas those greater than 0.5 and less than or equal to 1 represent non-swirling-like or 
hyperbolic regions. The value of 0.5 represents then a transition (parabolic) region where the magnitude 
of rotation rate and deformation rate are alike. The only exception applies to the Xcr/^d- Xcr/Xd- 
criteria. Since these criteria apply for elliptical regions only, we decided to normalize them from 0 to 0.5, 
where 0.5 is the boundary of elliptical region. 

The difference from the non-objective and objective versions lies on the reference frame from which the 
rate-of-rotation is computed. While the non-objective quantities use the original fixed frame, the objective 
quantities have their reference on the local frame attached to the eigenvectors of D. 

Normalized criteria are identified by a superscript star and are given by the following equations [13] 


Q* = — cos 

TT 


-1 




A* = — cos ^ 

TT 


\WP-\\D 


WP + \\DP 

2 \ 3 

' ^ 


fdet{D + W) 


21 


\W\\ 


\\D 


12 /det{D + W)y 


+ 


V 


V J 


AJ = 1-cos 


1 

TT 


-1 


Ao 


tr{D^) - Ai 




Al¬ 


lan 


-1 


A. 


DdW 


n 2 


D+W 




A* = — cos 

TT 


-1 


TT 

\wP 


WDP 


\D\\ 


^ fdet{D + W)d'' 


\WP + \\D\\ 


An = 1 -COS ^ 

TT 


^ fdetiD + W)d 


A 




( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 
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Al- 


tan 


-1 \ 


D+W 


TTT\ n 2 


D^W 


(18) 


3.1. Transient ABC flow 

The classic ABC flow [14,15,16] is largely used in the study of chaotic trajectories. Aiming to investigate 
high-frequency instabilities [17,18], we used a transient version of the ABC flow, also considered by Haller 
[2]. The flow field is given by the following set of equations: 


u{y, z, t) = A{t) sin( 2 ;) + B cos{y ), 

u(x, zfl)=B sin(x) + A{t) cos( 2 :), (19) 

w{x^ yfl) = C sin(^) + B cos(x). 

In the set of equations (19), A is time-dependent and is defined as A{t) = Aq + (1 + e~^^) sin(cc;t) with 
Aq = q = OA and uj = 27r, while B = \/2 and C = 1. 

Figure 1 contains the normalized classic flow classification criteria (Q*, A*, A 2 and on the first 

line and their respective normalized objective version (Q*, A*, A 2 and A*^/A*J on the second line. The 
non-objective criteria were obtained using the instantaneous velocity flow field described 
by the set of equations (19), whereas the objective criteria need two consecutive velocity 
fields to compute the rate-of-rotation of the eigenvectors of iA, (see Eq. (6) above). 



Figure 1. Iso-contours of the normalized flow classification criteria applied to the ABC flow fleld: non-objective (first line) 
and objective (second line) versions of Q*-criterion, A*-criterion, A^-criterion and /A*-criterion. 

The figure shows three surfaces of the cube whose dimensions are limited to the interval [0,27r]. Gen¬ 
erally, the normalized Q- and A 2 -criteria (first and third columns in Fig. 1, respectively) present very 
similar results, which may be explained by the strong relation between these two criteria (see [3]). The 
classifications of the normalized A-criteria (second column in Fig. 1) are also qualitatively similar when 
compared to those obtained using the normalized Q-criteria, although slightly different quantitatively, 
which, once again, may be justified by the relation between the A-criterion and the Q-criterion ([3]). Both 
A*^/A*- and A*^/A*-criteria (last column in Fig. 1) have a different range due to less detailed information 
on non-swirling-like regions. Nevertheless, they maintain a similar qualitative behaviour in terms of the 
location of vortex cores. 

Analysing the non-objective criteria (first line in Fig. 1) at the plane 2 : = 27r, the Q*-, A*- and A 2 - 
criteria (Figs, la-c) identify two vortex cores at {x ~ ~ '7r/2) and {x ~ ~ 37r/2), which are 
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characterized by the colour blue. Once again, due to its different range, almost all the region between 
0 < X < TT is blue according to the A*^/A*--criterion (Fig. Id), characterizing a swirling-like region. Still 
at this plane, all classic criteria seem to identify extensional (red) regions whose cores are located at 
approximately {x = 37r/2,^ = 7r/2) and {x = Znj2^y = 37r/2). The regions between elliptical (blue) and 
hyperbolic (red) regions are transition regions (green) generally referred to as parabolic^ where the role 
played by rotation and extension is equivalent. Qualitatively, the results at the plane y = 27r are very 
similar to those of the plane z = 27r. Two vortex (blue) cores are now identified at {x ~ 7r/2,z ~ 7r/2) and 
{x ^ 37r/2, z 7r/2), and two extensional (red) regions are located at approximately {x = 7r/2, z = 37r/2) 
and {x = 37r/2, z = 37r/2). It is worth noting that the magnitude of both swirling- and non-swirling-like 
regions are slightly greater than those identified at the plane z = 27r. Finally, for the plane x = 27r, 
the behaviour is yet similar, with vortex cores now identified at {y ~ 7r/2,z ~ 7r/2) and {y ~ 7r/2,z ~ 
37r/2), and extensional regions centred around {y ~ 37r/2, z ~ 7r/2) and {y ~ 37r/2, z ~ 37r/2). The main 
difference at this plane is the magnitude of the identified motions, which are reasonably weaker, indicating 
that the intensities of rotational and extensional motions are close to each other. 

From the perspective of the objective criteria (second line in Fig. 1), the flow at the three planes 
analysed above presents quite similar characteristics. The main difference is the remarkable increase in 
swirling-like (blue) regions. This fact is related to regions where the rate of rotation of the eigenvectors of 
D are more pronounced, see Fig. 2. This means that there are regions where the filaments aligned with 
the eigenvectors of D do not rotate intensively with respect to the reference frame where the problem is 
being described, but rotate significantly with respect to the frame attached to the eigendirections of D. 
This fact will be more explored in the following section for the 4:1 abrupt contraction below. 


3.2. Abrupt 4A contraction 

The flow generated by an abrupt 4:1 contraction presents different regions in which the fluid is submitted 
to shear, extension or rigid-body motion. Because this feature can be considerably useful for a rich 
discussion about flow classification, we considered a laminar steady state flow of a Newtonian fluid 
(Reynolds number based on the outlet velocity and height equal to 0.043). The numerical approach used 
is the same described by Mompean et al. [19] and the mesh consists of 150x80 grid points respectively 
in the streamwise (x) and wall-normal {y) directions. 

The iso-contours of non-objective (left column) and objective (right column) flow classification criteria 
are shown in Fig. 3. Likewise the ABC flow, the non-objective criteria employed here provide similar 
identifications and classifications except for the A*^/A*-criterion, since it has a different range associated 
to elliptical and parabolic domains only. The same rationale applies to the objective versions. 

Again, the green color represents regions which are not elliptical nor hyperbolic. In the present case they 
are related to the shear motion typical of fully-developed flow in a geometry of constant cross sectional 
area. 

The non-objective criteria Q*, A* and A 2 (Figs. 3a-c) identify a shear region away from the contraction, 
where the flow is fully developed. Close to the contraction and around the symmetry plane, an extensional 
region can be seen due the acceleration of the fluid passing through the contraction. Both near the sharp 
corner {x ~ 2.75,0.75 ^ ^ ^ 1) and the corner vortex {x ~ 2.5,0.1 ~ ^ ~ 0.25) a mix of extensional and 
rotational motions is observed. In fact, in these regions, the fluid is submitted to both extensions and 
rotations caused by the singularity of the geometry nearby. It is worth remembering that the A*^/A*- 
criterion does not provide any distinction in terms of non-swirling-like motions, i.e. swirling-like motions 
have a degree of intensity (greater than or equal to 0 and less than 0.5) whereas non-swirling-like motions 
possess only the classification (equal to 0.5, without a gradation). Therefore, in Fig. 3d, every region in 
red may be interpreted simply as a non-swirling-like region (possibly being either shear or extensional). 
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Figure 2. Eigendirections of tensor D in the three planes for the ABC flow considered. The ordering corresponds to the 
eigenvalues > X^ > X^. 

The blue regions, on the other hand, can be considered as elliptical regions. Thus, except for the size and 
intensity, the same two regions with rotational motions pointed by the other non-objective criteria are 
identified by the A*^/A*^-criterion. However, it also identifies two extra swirling regions: (a) the one close 
to the wall in the inner region, and (b) the one just by the outlet region. 

Concerning the objective criteria (Figs. 3e-h), the same three elliptical (swirling-like) regions identified 
by Mompean et al. [19] are observed here. As commented by the authors, the regions (i) just after the 
sharp corner {x 2.75,0.75 ^ ^ ^ 1), and (ii) near the corner vortex {x ^ 2.5,0.1 ~ ^ ~ 0.25) are 
intuitively understandable, since the flow is submitted to rotate at rates which are larger than the local 
rate-of-deformation. On the other hand, the third region (iii) just before the extensional region may be 
counter-intuitive at first, since it appears between regions of expected shear and extension. Again, the 
authors explain that this region should actually be considered as a plug flow (close to rigid-body motion) 
where the rates-of-rotation of the eigenvectors of D are larger than the local rates-of-deformation. One 
important result shown by the objective quantities is the information with a broad-spectrum of values for 
all the four criteria when compared with the non-objective quantities. This information is frame-invariant 
and related to the principal directions of the rate-of-deformation tensor, and is physically expected. 
Before entering a region where extension is dominant, the region just before the contraction plane, the 
eigenvectors of D rotate from the 7r/4 — 37r/4 directions of the shear flow in the fully-developed region to 
the 0 —7r/2 directions of the extensional flow. This rotation is overlooked by the non-objective criteria, but 
are captured by the objective ones. Although the filaments do not rotate predominantly with respect to 
the reference frame (x, ^), they do rotate with respect to the eigenvectors of iA, as clearly shown in Fig. 4. 
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Figure 3. Iso-contours of flow classiflers for the flow trough a 4:1 contraction: non-objective (left-hand side) and objective 
(right-hand side) versions of Q*-criterion, A*-criterion, A^-criterion and A*^/A*^-criterion. 


Because of this relative rate-of-rotation, new filaments are exposed to the corresponding eigenvalues. A 
similar rationale explains the appearance of the elliptical region just after the contraction, not detected by 
the non-objective quantities. The ^-eigenvectors need to rearrange in order to return to a fully-developed 
condition. 




Figure 4. Eigendirections of tensor D in the x — y plane for the 4:1 contraction flow. Again, the ordering corresponds to the 
eigenvalues > X^. 


It seems that the Acr/Aci-criterion has some problems on the so-called parabolic regions^ i.e. the regions 
of transition between hyperbolic and elliptical regions. It is difficult to delineate this intermediate region 
by this criterion in both, objective and non-objective versions. We can notice an almost binary color result. 
What we can see is that parts of the domain with no apparent physical difference that are classified as 
hyperbolic or elliptical. Most of these parts are in the parabolic region accordingly to other criteria. 
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4. Final remarks 


In the present work we analysed the performance of normalized objective versions of classic flow clas¬ 
sification criteria. The classic criteria are the Q-criterion proposed by Hunt et al. [4], the A-criterion 
proposed by Chong et al [5], the A 2 -criterion proposed by Jeong and Hussain [3], and the Xcr/Ki- 
criterion proposed by Chakraborty et al. [6]. The two flows considered were the transient ABC flow with 
A(t) = Ao + (1 + e~^^) sm{ujt) (with Aq = a/3, q = 0.1 and uj = 27r), B = a/2 and C = 1, and a 4:1 
abrupt contraction. The ABC flow is known for its chaotic character for the values employed. The abrupt 
contraction is a complex flow that exhibits different types of motion distributed in the domain: extension, 
pure shear, and rigid-body motion. 

The objective quantities in the present work use the effective vorticity as the entity associated to the 
elliptical character of the flow, in the place of the vorticity. The main difference between objective and 
non-objective quantities was the presence of elliptical regions in the objective versions which were not 
present in its non-objective counterparts. This happened because in these regions, the rate-of-rotation of 
the eigenvectors of the rate-of-strain tensor, with respect to the original frame, was significant and the 
effective vorticity uses this rate-of-rotation as reference for computing the rate-of-rotation of the filaments. 

Finally, it was shown that the objective criteria provides more information about the kinematics of the 
flow. Such new features of the flow may be useful in the investigation of complex flows and phenomena, as, 
for instance, drag-reducing flows, convection-driven problems or the mixture between two or more fluids. 
We emphasize that our goal is not to appoint to a preferable particular flow classification 
criterion, but to stand up for the advantages of applying criteria which enjoy objectivity. 
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